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A time domain subspace iteration technique is presented to compute a set of normal 
modes from the measured complex modes. By using the proposed method , a large 
number of physical coordinates are reduced to a smaller number of modal or 
principal coordinates. Subspace free decay time responses are computed using 
properly scaled complex modal vectors. Companion matrix for the general case of 
nonproportional damping is then derived in the selected vector subspace. Subspace 
normal modes are obtained through eigenvalue solution of the [M N ]~ l [K N ] matrix 
and transformed back to the physical coordinates to get a set of normal modes. A 
numerical example is presented to demonstrate the outlined theory. 


INTRODUCTION 


Experimental modal analyses are carried out to extract a set of modal parameters from the measured time or 
frequency domain data of the structure under test. These identified eigenvectors are in general damped complex 
modes due to several possible reasons t 1,2 !: 

1. The damping is nonproportional, i.e., [C] matrix is not proportional to [K] and [M] matrices. 

2. Measurement errors due to mass loading effects, noise, nonlinearities, etc. 

3. Digital signal processing errors due to finite frequency resolution, leakages, high modal density, and frequency 
response functions estimation procedure 

4. Modal parameter estimation errors due to invalid estimation of number of degrees of freedom. 

Identified complex modes can be used directly in the applications of modal modeling, structural dynamic modification 
and sensitivity analysis! 3,4 !, or validation and optimization of an analytical model! 5 ’ 6 !. On the other hand, real normal 
modes are sometimes more desired in the similar applications due to the facts that (i) normal modes are numerically 
easier to handle than complex modes and (ii) analysts usually computes normal modes rather than complex modes in 
the finite element analysis due to the lack of information of physical damping matrix. If normal modes are desired 
from a set of identified complex modes, a real-normalization procedure is needed to be implemented. Several methods 
have been proposed! 7 " 13 ^ the past to derive normal modes from identified complex modes. In this paper, a time 
domain subspace iteration technique (an improved approach from the method described in Reference 12) developed 
by the Structural Dynamics Research Laboratory (SDRL) is proposed to real-normalize identified complex modes 
originated from systems with nonproportional damping. 
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THEORY AND FORMULATION: in the physical coordinates 


Assume there exists a set of measured modal parameters from a modal test. These modal parameters consists 

of N complex modes (and their complex conjugates), A r , r=l, N and {^} r , r=l, •••,N. Each modal vector has 

dimensions N m and in general N m > >N. To compute a set of normal modes from the given complex modes, a time- 
domain approach developed in the physical coordinates is formulated. 


Formulation of Free Decay Time Responses 

For the given modal parameters, displacement, velocity and acceleration responses can be expressed as 


{*(0}= 

(1) 

r=l 

(2) 

r=l 

{i'(0}= %m}re Kt 

r=l 

(3) 


Computation of Companion Matrix [E] 


The equations of motion for the general case of nonproportional damping is 


[C] [M]iri(r)| 

m [0]Jji‘(r)j + 


\K] [onrxcoi r[fn 

[0] -[M]jj i( ,)} = j[0]} 


Another form of the homogeneous Equation 4 can be written as 



■ [0] [/] 1 
-[MHA] -[MHC] 


[■*(0 

m 


Let 


[£] = 


[ 0 ] [/] • 

-wnK] -wnc] 


(4) 


(5) 


( 6 ) 


Equation 5 can then be written as 


{y(0} = [£]{>-(0} 


(7) 


where {y(t)} is now the system’s state vector containing the displacement response (x(t)} and velocity responses (i(f)} 
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and [E] is called the companion matrix. In the physical coordinates, by repeating Equation (7) 2N m time instants and 
dropping vector notation, the following equation is satisfied 


tv(^i), • • • = [ £ ]Ly('iX • • • >y('2ArJ] 


( 8 ) 


where t x ,/ 2 , ■ ■ ■ ,t 2Nm represent 2N m time instants. From Equation 8, the companion matrix [E] can be computed 


[£]=m [Y ]- 1 


( 9 ) 


where [E], [Y] and [Y] are all 2N m x2N m matrices. 


Computation of Normal Modes 

By computing the companion matrix [E], [M] _1 [X] matrix gives the normal modes according to the eigenvalue 
solution 


[[MW]] [<*>]= w r 2 [<H r = 1, • • • ,N 


In the physical coordinate, [Y] is always singular due to the number of measurement degrees of freedom N m is much 
larger than the number of measured modes N. In other words, it is numerically difficult and unstable to soTve rEl in 

the physical coordina t e s • . Therefore, it is proposed to compute the companion matrix [E] in the principal or 
modal coordinates in the following way. 


NORMAL MODES SOLUTION: subspace iteration technique 


Based on a set of identified complex modes, a time domain subspace iteration technique is developed to obtain 
a set of normal modes in the following way: 


Hypothesis 


A 

Define [$] = transformation matrix of the identified complex modal matrix [$], it is assumed that f^l can be 
expressed as 1 J 


[*] = m [W] 


(ID 


where [W] is a NxN complex matrix and can be obtained through pseudo inverse technique 


m =[*]+[*] 


( 12 ) 
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For the proposed method, the selection of the transformation matrix [<f>] can not be arbitrary since it represents a N- 
dimensional subspace of the physical A^-dimensional vector space from which the companion matrix [E N ] and the 
normal modes are derived. Therefore, the following method to determine the transformation matrix [$] is proposed 
and justified. 


Determination of the Transformation Matrix [<f>] 


Using the proposed method, it can be shown that choice of the [<S] has significant effects on the calculated 
normal modes. Since the transformation matrix [$] can be considered as a N-dimensional vector subspace (column 
space of [$]), into which the physical system matrices and the identified complex modes are transformed, therefore, it 
must be so chosen such that the following conditions are satisfied: 

1. [$] must not be orthogonal to the normal modal matrix [$] of the undamped system, and 

2. [&] should be selected that it is as close to [<J>] as possible. 

3. [$] must has rank N. 

A logical way to select [$] is based on the identified complex modal matrix ['£]. 

1. Fillod^ 14 ^ proved that if the complex modal vector is normalized according to the following equation 

2 xmrwmr + wncm = 2 m a3) 

,the imaginary part of {ip} r i s minimized. This indicates that the real part of {rp} r is maximized and contains 
maximum useful information of the identified complex modes. 

2. It can be proved (see Appendix) that if a undamped system is perturbated to the first order by a 
nonproportional damping matrix, and the calculated complex modes are normalized according to Equation 13, 
then the real parts of the normalized complex modes are very close to the normal modes of the undamped 
system. 

3. The transformation matrix obtained by taking the real part of a set of complex modes normalized according to 
Equation 13, has rank N due to the fact that (i) the identified complex modes are independent of each other, 
and (ii) norms of the columns of [l>] have the same order of magnitude, i.e, eigenvalues found in Equation 19 
using singular value decomposition technique also have the same order of magnitude. 

From the above observations, it can be concluded that real part of the complex modes normalized to 2 ju r according to 
Equation 13 is the best choice of the transformation matrix [$]. 


Computation of Free Decay Time Responses in Subspace 

In order to derive the companion matrix [£*] in the selected subspace [4>], free decay time responses in the 
physical coordinates needs to be transformed into the selected vector subspace using one of the following two 
approaches: 

Modal Space Approach 

The displacement vectors in the physical and modal coordinate systems are related by the following 
transformation: 
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{*(/)} = [$]{/>(/)} 


( 14 ) 


or 

{pit)} = t$r woj (is) 

where [$] + is the pseudo inverse of [$]. 

Substitute Equation 1 and 11 into 15, for the given modal parameters, free decay displacement in the modal 
coordinate can be written as 

iP(t)}= [*r2Re{[®]{e M }} 

= 2 [$]+[*] Re{ [IT] {* M }} ( 16 ) 

= 2Re{[»']{e v }} 

where "Re" represents the real part of a matrix. 

Similarly, velocity and acceleration responses in the modal coordinates can be expressed as 


{/>(')}= 2Re{[^]{Ae v }} 

(17) 

[Pit)} = 2ReUW']{A 2 e A ' t }} 

(18) 


Principal Response Analysis (PRA) Approach 


Using singular value decomposition technique and assuming [$] has full rank N, [<£] can be decomposed as 


[*] = [P] f Ej [Sf 


(19) 


where 


[P] = orthonormal matrix (N m xN) 
r Ej = diagonal matrix consists of 

eigenvalues of [$]*[$] (NxN) 

[S ] = unitary matrix consists of 

eigenvectors of [$]*[<£>] (NxN) 

Using matrix [P] as the transformation matrix, displacement in the physical coordinates can be transformed into the 
principal coordinates according to 


«')} = [P]{R(0} 


( 20 ) 





(21) 


Substitute Equation 1 and 19 into Equation 21,{p(t)}can be expressed as 


{p(t)}= [P] ff 2ReUtf]{e M }} 

= 2[?f[4]Re{[^]{e M }} 

= 2[Rf[P]rE J [SfRe{[^]{e Art }} 

= 2 rSjtSfReUW'Ke^}} ( 22 > 

Similarly, velocity and acceleration responses can be written as 

(p(0> = 2 fSj [SfReUlPKAe^}} (23) 


{p(t)} = 2 TEj [S] 1, Re{[^]{A 2 e Art }} (24) 


Computation of companion matrix [ E N ] 


Similar to the derivations (Equation 4 to Equation 6) in the physical coordinates, the companion matrix E N can 
be computed in the selected vector subspace using the time responses in equations 16-18 (modal space approach) or 
22-24 (PRA approach) 


{<7(0} = [£*](<)} ( 25 > 

where (q(t)} is now the system’s state vector in the selected subspace containing the displacement response (p(t)} and 
velocity responses {p (7)}, and 


[E*] = 


[ 0 ] [/] ‘ 


(26) 


where [M N ], [X N ] and [C*] are the reduced system mass, stiffiiess and damping matrix respectively and can be written 
as 


[MjrW = [if [M] [$] (modal space) 

(27a) 

[M w W = [Pf[M][P] (PRA) 

(27b) 

= (modal space) 

(28a) 

ffirW-Pl'WP] ( PRA ) 

(28b) 

[C w W = [4>f[C][i] (modal space) 

(29a) 

[^W = [Pf[C][P] (PRA) 

(29b) 
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Once {q (/)} and {4(0} are computed at 2N time instants, companion matrix [E n ] Wx2 n can then be computed in the 
modal or principal coordinates accordingly 


[4(0)» ‘ ‘ ‘ »4(Oj0]-[£w]14(O)> ' • ’ >9(Ow)] 


(30) 


and 


[£*]= mm 1 


(31) 


Since [Q ] is a matrix with full rank 2N, it is always invertible. 


Computation of Normal modes: first iteration 

From the companion [£y], the N by N [AT] 1 ^] matrix gives the normal mode solutions {$] in the modal or 
principal coordinates according to the eigenvalue solution 


[[M„m„]] © = a 2 ffl, r = 1, ,N 


(32) 


where ft r is the r-th undamped natural frequency. 

From the modal space method, it is noticed that (i) if [4] is selected as the normal modal matrix [$] of the undamped 
system, then [Afy], [£y] and consequently, [Afy]* 1 ^] are all diagonal matrices by the orthogonality conditions, and 
(ii) Generalized damping matrix [C^ n ] can also be obtained from the companion matrix [£y] as [My ]’ 1 [Ay], 
Generalized damping matrix is defined as 


{ c ^] = m T [cjw 

and [$] must also satisfies the following relationship 

m T m m = [/] 


(33a) 


(33b) 


First estimation of the normal modes {^>} r ,r=l, • • • ,N of the undamped system can then be obtained through the 
coordinate transformation 


[*ii= mm (modal space) 
or 

[$i] =[£][£] (PRA) 

where is the first estimation of the normal modal matrix. 


(34a) 


(34b) 
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Subspace Iteration Technique 


In order to improve (if possible) and check the accuracy of [$i], [fy] is compared with [$] column by column 
to check the convergence. If not convergent, then [$i] will substitute [$] in Equation 11,14 and 19 and a new iteration 
will start to find a [$ 2 ]- This process continues until [$ m ] converges to where m is the number of iterations. In 

general, m is usually 2 or 3, i.e., the transformation matrix [$] converges very fast because all iterations are done in the 
same vector subspace which is defined by the initially estimated transformation matrix [$]. 


Y 


0 



Figure 1. Folded Beam Used as Test Structure 


CASE STUDY 

The proposed method is applied to the identified complex modes of a simple structure as shown in Figure 1 
(For comparisons, this example is the same as the one presented in Reference 11). This structure is a folded beam 
excited on its bending modes along the y axis. Its dynamic characteristics shows the existence of pairs of eigenvectors 
at very closely spaced eigenvalues (quasi multiplicity of order 2). The eigenvectors are described by seven 
measurement degrees of freedom as indicated in Figure 1. The first two complex modes are identified from a set of 
frequency response functions measured using slow sine sweep excitation technique. The identified complex modes 
scaled according to Equation 13 are listed in Table 1. The companion matrices [E N ] calculated in the selected 
subspace using both the modal space and the PRA approaches are listed in Table 2. Two undamped natural 
frequencies and their corresponding normal modes based on the proposed method and the three methods proposed in 
Reference 11 are listed in Table 3. Modal Assurance Criteria (MAC) [15] are computed for these two normal modes 
between those obtained from the proposed method and those obtained from the methods described in Reference 11. 

From Table 3, both the modal space and the PRA approaches calculate the same undamped natural frequencies and 
normal modes of the associated undamped system. 

From Table 2, companion matrix [E N ] obtained from the second iteration of the modal space approach shows (i) the 
Wn]' 1 [Kw] matrix is fairly diagonalized and (ii) the [Afy]' 1 [C N ] matrix is fairly symmetrical. These indicate that the 
calculated normal modes are very close to the true real modes of the associated undamped system. Therefore, the 
calculated [Af^]' 1 [C w ] matrix listed in Table 5 is also very close to the generalized damping matrix derived in 
Reference 11. 

From Table 3 and 4, two normal modes derived from the proposed method stay very close to those derived in 
Reference 11. 
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TABLE 1. Measured Complex modes 


Measured Complex modes 


First Mode 

Second Mode 

Point No. 

Ai = -1.992 + j 209.55 

A 2 - - 

2.547 + j 213.40 


real 

imaginary 

real 

imaginary 

1 

-0.144 

0.118 

0.266 

.0683 

2 

-0.303 

0.203 

0.457 

0.138 

3 

-0.498 

0.267 

0.598 

0.228 

4 

-0.179 

0.231 

0.519 

0.0847 

5 

0.210 

0.207 

0.464 

-0.0942 

6 

0.618 

0.192 

0.433 

-0.280 

7 

1.060 

0.192 

0.416 

-0.483 


TABLE 2. Calculated Companion Matrix ] 


Companion matrix [£#] 

Iteration No. 

Modal Space Approach 

PRA Approach 

1 

0.0000 0.0000 1.0000 0.0000 
0.0000 0.0001 0.0000 1.0000 
-43507.2 -120.312 -3.7156 4.3521 
-117.281 -45973.5 4.3826 -5.3629 

0.0000 -.0001 1.0000 0.0000 
0.0000 0.0000 0.0000 1.0000 
-43970.09 -1182.15 -0.8233 2.9472 
-796.414 -45510.59 2.0166 -8.2553 

2 

0.0000 -.0001 ' 1.0000 -.0000 

0.0000 0.0000 0.0000 1.0000 

-43501.49 0.0000 -4.1383 4.3492 

-0.0156 -45979.18 4.5045 -4.9402 

-.0000 0.0001 1.0000 0.0000 
0.0000 -.0000 0.0000 1.0000 
-43921.51 -1142.37 -0.9521 3.1301 
-756.586 -45559.18 2.1995 -8.1294 


TABLE 3. Calculated Normal Modes 


Calculated Normal Modes 


First Mode 

Second Mode 

Method 

1 

Method 

2 

Method 

3 

Modal 

Spa. 

PRA 

Method 

1 

Method 

2 

Method 

3 

Modal 

Spa. 

PRA 


Tic 

itural fre 

quency ( 

rad/sec; 


Tic 

itural fre 

quency ( 

rad/sec; 


Point No. 

208.38 

208.59 

208.68 

208.58 

208.58 

214.22 

214.45 

214.32 

214.43 

214.43 

1 

.161 

.153 

.153 

.151 

.151 

.458 

.454 

.454 

.451 

.451 

2 

.333 

.312 

.318 

.312 

.312 

.774 

.769 

.771 

.771 

.771 

3 

.535 

.506 

.515 

.506 

.506 

1.0 

1.0 

1.0 

1.0 

1.0 

4 

.215 

.198 

.200 

.196 

.196 

.905 

.890 

.892 

.889 

.889 

5 

.169 

.179 

.178 

.181 

.181 

.867 

.830 

.836 

.827 

.827 

6 

.571 

.572 

.577 

.574 

.574 

.873 

.810 

.826 

.815 

.815 

7 

1.0 

1.0 

1.0 

1.0 

1.0 

.901 

.82 

.832 

.815 

.815 
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TABLE 4. MAC for the Calculated Normal Modes 


Modal Assurance Criteria 


First Mode 

Second Mode 

Method 1 

Method 2 

Method 3 

Method 1 


Method 3 

PRA 

.999122 

.9999906 

.9999941 

.998605 

MM 

.999913 


TABLE 5. Generalized Damping Matrix [C^] 


Method 2 

Method 3 

Proposed Method 

4.86 -4.72 

-4.72 5.50 

4.14 -4.27 

-4.27 4.94 

4.14 -4.35 

-4.50 4.94 


CONCLUSION 

The proposed method described in this paper has the following advantages: 

1. This method is numerically very efficient and stable because all computations are performed in the selected 
vector subspace and it requires very few iterations. 

2. Since this method uses a time domain approach, it can be proved that the calculated normal modes are not 
sensitive to the norm or scaling errors existing in the identified complex modes. 

The drawbacks of this method are: 

1. It can not improve the norm or scaling errors existing in the identified complex modes. 

2. The calculated normal modes are subject to modal truncation errors since an incomplete set of modes are 
identified. 

It can be concluded that the success of the proposed method is dependent on the selection of the transformation 
matrix [$] and the quality of the measured complex modes. It seems to be a reliable method that can be used as a 
post -processing procedure of the identified complex modes for further applications. 
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APPENDIX 


The characteristics equations of a undamped and a damped system can be written as 

[*][*] -[M][*] W j-Q. (35) 

[M][tt] r A 2 j + [C][¥] fAj + [*][*] = 0 (36) 

It can be proved that in the case of a lightly damped system, i.e, [C] = e [C], where e is a small perturbation parameter, 
the natural frequencies and the normal modes of the associated undamped system can be approximated by taking the 
damped natural frequencies and the real part of the complex modal vectors respectively. 

Assuming: 


[*F[cp] = [r] 

(37) 

r Aj = j faj + e If/, 

(38) 

M = [$][[/] +e ,[Fj 

(39) 
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and recalling that complex modal vectors satisfy the following orthogonality conditions: (Note that both the Equation 
13 and 40 are scaled to 2 /' w) 

[*F[M][tt] r Wj + r Wj [®f [Af ][®] + e[^] T [C][^] = 2 j fwj 
[*‘FM*] r «j + r A*j[^F[M][^] + «[*T[C]M = o 

where represents the complex conjugate of a matrix. 

Substitute Equation 35,36 and 37 into Equation 40 and 41, equating the first order terms of e: 


[r] = -2 rt/j + faj [v + v T ] + [v + e 7 "] rnj (42) 

[r] = -2 rt/j + iQj [k - v T \ + [k - f 7- ] rnj (43) 

and then 

[r] = -2 rt/j + [F] faj + r nj[FF ( 44 ) 

o= rftj[F] + [FF r fij ( 45 ) 


From Equation 45, the matrix r flj[F] is anti-symmetrical, therefore the diagonal elements of this matrix are all zero. 

Since [r] is a real matrix. From Equation 44, it can be concluded that: 

1. f[7j is a real diagonal matrix. This indicates that the effect of the first order perturbation damping matrix on the 
eigenvalues is purely real. 

2. [V] is a real matrix. This indicates that the effect of the first order perturbation damping matrix on the 
eigenvectors is purely imaginary. 


(40) 

( 41 ) 
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